Voltage-driven quantum oscillations in graphene 
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We predict unusual (for non-relativistic quantum mechanics) electron states in graphene, which 
are localized within a finite-width potential barrier. The density of localized states in the sufficiently 
high and/or wide graphene barrier exhibits a number of singularities at certain values of the energy. 
Such singularities provide quantum oscillations of both the transport (e.g., conductivity) and ther- 
modynamic properties of graphene — when increasing the barrier height and/or width, similarly 
to the well-known Shubnikov-de-Haas (SdH) oscillations of conductivity in pure metals. However, 
here the SdH-like oscillations are driven by an electric field instead of the usual magnetically-driven 
SdH-oscillations. 

PACS numbers: 73.22.-f, 



The Shubnikov-de-Haas effect, i.e., the oscillations of 
the magneto-resistance of metals when increasing an ex- 
ternal magnetic field, was one of the first macroscopic 
manifestations of the quantum-mechanical nature of mat- 
ter. The key to understanding this remarkable phe- 
nomenon was pointed out by Landau and Onsager and 
it is described in many textbooks on solid state physics 
(see, e.g., Ref. [l]). Namely, electrons in the conduction 
band of a metal in a strong magnetic field behave like 
simple harmonic oscillators. The resulting energy spec- 
trum is made up of equidistant Landau levels separated 
by the cyclotron energy. The density of electron states 
has singularities at the Landau levels. When the mag- 
netic field is changed, the positions of the Landau levels 
move and pass periodically through the Fermi energy. 
As a result, the population of electrons at the Fermi sur- 
face also oscillates and, in turn, leads to quantum oscilla- 
tions of the transport and thermodynamic properties of 
a metal. The quantum oscillations also manifest them- 
selves in the thcrmoconductivity, magnetization, sound 
attenuation, magnetostriction, and other quantities. 

These quantum oscillations are pronounced in con- 
ductors with a long mean free path of charge carriers. 
This can occur in pure metals, scmimctals, and narrow 
band-gap semiconductors at low temperatures, as well 
as in graphene, a one-atom-thick sheet of carbon. The 
Shubnikov-de-Haas oscillations of the magneto-resistivity 
were observed in graphene 0, Q soon after its discov- 
ery Due to the monolayer honeycomb-lattice struc- 



ture of graphene, its electrons obey a masslcss Dirac-likc 
equation (see, e.g., Refs. H, @). This is responsible for 
the unusual properties of graphene. In particular, the 
Landau levels in graphene are not equidistant and these 
influence the period of the Shubnikov-de-Haas oscilla- 
tions @, S|- Graphene has another striking property: it 
has unusual relativistic effects which are counterintuitive 
for electrons with speeds much slower than the speed of 
light For example, it has been recently shown [5| that 
graphene could be used for experimentally testing the 
so-called Klein paradox [Io| . This quantum- mechanical 
effect of relativistic particles penetrating through high 
and wide potential barriers can be illustrated with mass- 
less Dirac fermions in graphene with a potential barrier 
controlled by an applied voltage. A high pot ential energy 
barrier in graphene, as was shown in Ref. [ll|, can also act 
as an unusual electron lens, due to the negative refraction 
of electron waves at the edge of the barrier, in analogy 
to the negative refraction of 3D [l^] and 2D [13 , 14 1 elec- 
tromagnetic waves. 

Our goal here is to show that, due to the Dirac-like 
Hamiltonian of graphene with a potential energy barrier, 
quantum oscillations similar to the Shubnikov-de-Haas 
effect can be observed without an applied magnetic field. 
Below we prove that the density of electron states in a 
graphene sheet with a potential barrier should display 
quantum oscillations if the strength of the barrier (i.e., 
the product VqD of the barrier's height Vb and width 
D) exceeds some threshold value. In these oscillations, 



the barrier strength VqD, that can be controlled by a gate 
voltage, plays the same role as the external magnetic field 
in the Shubnikov-de-Haas effect. 

The quantum oscillations predicted here originate from 
a new type of electron states in graphene. Contrary 
to non-relativistic quantum mechanics, where localized 
states can only exist inside quantum wells, we find that 
the electron states in graphene can be localized within 
the barrier. We show that the energy E(q y ) of the local- 
ized states (versus the wave vector component q y along 
the barrier) becomes non-monotonic if VqD > ttHvf 
[vf is the Fermi velocity). This produces singulari- 
ties of the density of electron states for energies where 
dE/dq y = 0. When the magnitude and/or width of the 
barrier changes, the locations of the singularities move 
and periodically cross the Fermi level, generating quan- 
tum oscillations of both thermodynamic and transport 
properties, e.g., of the conductance in the y-direction 
(along the barrier). 

Electron states localized in a barrier. — The tunnelling 
of relativistic particles through a finite-width potential 
barrier has recently been studied in Refs. HE IS Here 
we consider another type of electron waves that propa- 
gate strictly along the barrier and damp away from it. 
Our analysis shows that a step- like barrier (i.e., a sin- 
gle edge of an infinitely-wide barrier) does not support 
such electron waves, which would be an analog to sur- 
face electromagnetic waves (plasmon-polaritons) at the 
interface between two different media. Therefore, even 
though the potential barrier in graphene could act as an 
electron lens fill ], it cannot provide the perfect lensing 
(i.e., subwavelength image reconstruction) that is pos- 
sible for Veselago's lens in optics [l7j]- However, as we 
show in this section, the electron waves in graphene can 
be localized inside a finite-width potential barrier. 

We consider electron states in graphene with a poten- 
tial barrier located in a single-layer graphene occupying 
the xy-plane (see Fig. 1). For simplicity, we assume that 
the barrier V(x) has sharp edges, 

' 0, \x\ > D/2, 
V , \x\ < D/2. 

Electrons in monolayer graphene obey the Dirac-like 
equation, 
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Hip = ih 
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where vp is the Fermi velocity and a = (a x , cr y ) are Pauli 
matrices. We seek stationary spinor solutions of the form, 



ip = ip(x) exp 



-Et + iqy 



(3) 



with energy E and momentum hq along the barrier. We 
focus on the states with 




x/D 



FIG. 1: (Color online) (Top) Geometry of the problem. A 
graphene sheet is placed under the voltage gates indicated by 
block rectangles. (Bottom) Potential energy barrier V(x) in 
graphene (dashed blue line) and the probability distribution 
W(x) — \ipi(x)\ 2 + \ip2(x)\ 2 for the localized electron state at 
qd = 3.5625, eD = 0.003, and V = 9 (red solid line). 



In this case, the electron waves satisfying Eq. @ damp 
away from the barrier, and the components ip\ and "02 of 
the Dirac spinor can be written in the form, 



aexp(k x x), x<—D/2, 
bexp(iq x x) 

+cexjp(-iq x x), \x\ < D/2, 

dexp(—k x x), x > D/2, 



— b exp(iq x x + iff) 
+cexp(—iq x x — iff), \x\ < D/2, 
Tj ^| J exp(-fc x z), x>D/2 



(5) 



(6) 



with real k x = (q 2 -e 2 ) 1 / 2 and q x = [{e-V / D) 2 - q 2 ] 1 / 2 . 
Here V = VqD/Hvf is the effective barrier strength and 
tan = q/q x . 

Matching the functions ipi(x) and ip2(x) at the points 
x = ±D/2, we obtain a set of four homogeneous algebraic 
equations for the constants a, b, c, and d. Equating the 
determinant of this set to zero, we obtain a dispersion 
relation for the localizes electron states, 



teLii(q x D) 



k x q x 



(V/D - e)e + q 2 



(7) 



\q\ > \e\ = \E\/Hv F . 



(4) 



Figure 1 illustrates the behavior of the probability distri- 
bution W(x) for a localized state. Note that W(x) is an 
even function with continuous derivative W'{x), in spite 
of the fact that each one of the functions \tpi\ 2 and \ipi\ 2 
are not even (the chirality of the Dirac spinors) and both 
have discontinuous derivatives at the points x = ±D/2. 

Note that similar localized states can also be observed 
in a 2D electron gas when a voltage is applied to produce 
a potential well. For electrons with a quadratic dispersion 
law, this spectrum is £ n = y 2 + AQ 2 — V, where £ = 
mD 2 E/2h 2 , V = mD 2 V /2h 2 , Q = qD, m is the electron 
mass. Here y n is the nth root of the equation ytan(j/) = 



3 



(V — y 2 ) 1 / 2 for even states, and ycot(y) = — (V — y 2 ) 1 ^ 2 
for odd states. 

The spectrum Eq. ([7]) of localized states in graphcnc 
is shown by the solid black curves in Fig. 2, for dimcn- 
sionless variables Q = qD and £ = eD. This spectrum 
consists of an infinite number of branches £ n (Q). Each 
of these branches starts from the lines £ = i\Q\ (red 
solid straight lines in Fig. lb) at £ = V/2 - 7r 2 n 2 /2V 
and tends asymptotically to the lines £ = V ± Q (dashed 
red lines). Moreover, a particular branch of the spectrum 
starts at the point (Q = 0, £ = 0) and also tends to the 
line £ = V — Q, when increasing Q. 
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FIG. 2: (Color online) Electron spectrum in graphene ob- 
tained for V = 1 (inset) and V = 9 (main panel). The sea 
of delocalized states (continuum spectrum) is marked by the 
light-purple regions. The branches of the spectrum for lo- 
calized states are shown by solid black curves between the 
straight solid and dashed red lines. There are no states in the 
forbidden (white) regions. 

The behavior of different branches of the spectrum 
depends on the barrier strength V. If V < 7r/2, all 
branches satisfy £ < (see inset in Fig. 2). Localized 
states with positive energies appear only for V > n/2. 
When V increases, new branches in the spectrum with 
positive energies appear. When V is within the interval 
(??. + l/2)7r < V < (n + 3/2)ir, the number of branches 
with E > is (n + 1), n = 1, 2, 3, We empha- 

size that each of the branches with positive energy has 
a maximum £™ ax at a certain wave number Q = Q™ ax . 
Near these points, the group velocity of localized elec- 
tron waves tends to zero. This effect is similar to the 
stop-light phenomenon 11811 found in various media, in- 
cluding superconductors [19j . 

Note that defect-induced localized electron states in 
graphene and the enhancement of conductivity due to an 
increase of the electron density of states localized near 



the graphene edges were recently reported 2(], 21, 22[. 
Contrary to these examples, the electron states studied 
here are localized within the barrier and also these are 
tunable, i.e., the energy levels can be shifted by changing 
the barrier strength (e.g., via tuning a gate voltage). 

Density of localized electron states. — To calculate the 
density N(E) of electron states, we use the general for- 
mula N(E) — J2 a ~ E a ), where a labels the quan- 
tum state and 5(x) is Dirac's delta-function. Using 

T ■■■ 

spectrum, we derive 



2L x L v {2'k) 2 J_ oo dk x dk y . . . for continuum 



N \£\, N Q 



L T .Li 



(8) 



TrhvpD ' 

where L x and L y are the lengths of the graphene sheet in 
the x and y directions, respectively. For localized states, 
we obtain 

-l 



N loc (£)=2N - 
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d£ n (Q) 



dQ 



(9) 



e n (Q)=s 



Here n runs over the number of positive roots of the 
equation £{Q) = £. 

The dimensionless density of states N(£)/Nq is shown 
in Fig. 3. The localized electron states exhibit two types 
of peculiarities. First, increasing £, the jumps or steps 
(each one of magnitude 2D/L X ) in N(£)/No occur at the 
points £ = V/2 — ir 2 n 2 /2V, where new branches of the 
spectrum arise or disappear. More importantly, singu- 
larities are observed when £ = £™ ax , where {d^/dQl" 1 
in Eq. ((9]) diverges. 
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FIG. 3: (Color online) Dimensionless density of electron states 
N(£)/Nq in graphene with a potential barrier, for D / L x — 0.5 
and V = 16. 

The locations of the singularities shift when changing 
the barrier strength V. Therefore, they periodically cross 
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the Fermi level £p- This produces quantum oscillations 
of the density of states at the Fermi energy. They are 
seen in Fig. 4, showing N(£p)/N n versus the effective 
barrier strength V. 
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FIG. 4: Dimensionless density of electron states N(J\f)/No at 
the Fermi level versus the effective strength V of the potential 
barrier, for D/L x = 0.5 and £f = 1. 

The periodic change in the number of electron states 
near the Fermi level, increasing the barrier strength (e.g., 
by varying a gate voltage), necessarily results in quan- 
tum oscillations of the transport and thermodynamic 
properties of graphene. For example, the conductance 
of graphene along the barrier qualitatively mimics the 
quantum oscillations of the density of electron states. 

Figure 4 shows the quantum oscillations of the den- 
sity of electron states in e-typc graphene (with a positive 
Fermi energy). For p-type graphene, with Ep < 0, quan- 
tum oscillations of the density of states at the Fermi level 
can also be observed, if the (now opposite-bias) applied 
voltage forms a potential well instead of a barrier. Indeed, 
the Dirac equation ^ is invariant with respect to the 
transformation: E — > — E, V — ► —V, x —> —x, y — ► —y. 

In conclusion, we predict an unusual type of electron 
states in graphene localized within a potential barrier. 
For barriers with sufficiently high magnitude and width, 
the density of localized states has singularities. This fea- 
ture of localized states can result in quantum oscillations 
of the thermodynamic and transport properties (e.g., the 
conductance along the barrier) of graphene when chang- 
ing the barrier strength (e.g., by varying a gate voltage). 



These electric-field driven quantum oscillations are simi- 
lar to the Shubnikov-dc-Haas oscillations of conductivity, 
which are produced in standard metals when changing 
the external magnetic field. 
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